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ABSTRACT 

This work is a continuation of plasma oscillation studies 
made in rare gas plasmas extending over the last four years at the 
Naval Postgraduate School. An electronic method of obtaining con- 
tinuous graphical plots of probe current vs. applied probe potent- 
ial and probe current Squared vs. applied probe potential was 
developed and successfully used. Preliminary experiments indicate 
this system may be extended to obtain the first and second deriva- 
tives of probe current with respect to applied probe potential 
allowing an accurate determination of space potential and.electron 
energy distribution functions. A gas diode filled with argon at 
5.9 microns of mercury was used in this investigation. Probe data 
was taken for twenty positions along the axial length of the tube 
using an anode-cathode spacing of 42 mm. Five different discharge 
currents were studied. A second run was made for ten positions 
and a discharge current of 50 ma at 6.9 microns of mercury to study 
distribution functions. Electron oscillations were detected which 
corresponded with the Langmuir plasma frequency predicted from the 
plasma density; the appearance of these oscillations corresponded 
with a decreasing electron density gradient. Peak intensity of 
oscillations occurred along a decreasing density gradient. 
Oscillations were also measured in the kilocycle range which are 
attributed to ion oscillations associated with beam electrons. 
Data iS presented in the form of: raw data plots; space potential 
as a function of probe position; electron temperature as a function 


of probe position; electron density as a function of probe position; 


at 





electron density as a function of discharge current; and electron 
energy distribution functions. It is tentatively suggested that 
the unmodulated beam of charged particles, coming from the cathode, 
passing through and maintaining the plasma, suddenly becomes un- 
stable and gives rise to electron oscillations of the plasma in 


the presence of an associated decreasing plasma density gradient. 
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l. Introdwe@tien: 

The recent successes in harnessing the explosive fusion of 
the nuclei of the Light elements and the even more recent attempts 
at space exploration has led to a renewed peak of interest in 
plasma physics in an effort to better understand the basic physical 
processes occurring in a highly-ionized, gaseous, conducting medium. 

Success in understanding these phenomena may ultimately lead 
to generation of fusion power through the medium of a plasma. 
Linhart (1) describes this possibility in the following way. 

It was appreciated that the energy output from fusion 

reactions depends critically on the kinetic energy of 

the colliding nuclei and that fusion outputs of 

practical interest depend upon one's ability to produce 

temperatures of at least several million degrees Kelvin. 

If explosions are to be avoided, then the pressure of 

matter at this temperature must be balanced by external 

forces. This is within the power of our engineers only 

if the density of the nuclear fuel is less than the 

density of our atmosphere. The search for a mechanism 

of a controlled release of fusion energy became, there- 

fore, synonymous with the study of high temperature, 

low density plasmas. 

Since 99.9% of the matter in our universe is ionized and 
therefore composed of plasma, Communication and electronic systems 
used in space exploration must be able to operate in, and propagate 
Signals through, ionized regions of the outer atmosphere and space. 


Thus, an understanding of plasma processes will contribute strongly 


to the solution of these space problems. 
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Plasma physics is concerned with the behavior of systems of many 
free electrons and ionized atoms where the mutual Coulomb interactions 
cannot be ignored. Ina restricted sense, such systems of particles 
consist of nearly equal numbers of positive ions and negative 
electrons. Systems of this type are examples of a medium known as 
plasma which in many respects behaves differently from the solid, 
liquid and gaseous states of matter (1). 

The average kinetic energy per particle determines the state in 
which matter exists. Matter can exist as a plasma, 1.e., in its 
fourth state, when the kinetic energy per plasma particle exceeds 
the ionizing potential of atoms -- which is usually a few electron 
mikes (1). 

Plasma spans a broader energy band than any other state of 
matter; it occupies about 20 octaves of the kinetic energy spectrum. 
This wide kinetic energy band of the plasma state is the reason for 
so much common ground between plasma physics and many other fields of 
physics, such as the dynamics of single charged particles, or the 
physies of electrical discharges in gases (1). Linhart says: 

Although probably more than 99.9% of the matter in our universe 

is ionized and therefore in the plasma state, on our planet 

plasma has to be generated by special physical processes and 
special conditions. 

On earth the medium which most closely resembles an ideal plasma 
is a partially ionized gas. This fourth state of matter defined as a 
plasma is not very much at home on our dense, cold planet; therefore, 
Special conditions must be realized in order to generate plasma-Llike 


conditions in the laboratory. These special conditions are in the 


realm of the physics of electrical discharges in gases (gaseous 


electronics) and hence, the close kinship between this discipline and 


plasma physics. 





Plasma was originally defined by Langmuir and Tonks (2) 


in the following way: 

The word 'plasma' will be used to designate that portion 

of an arc type discharge in which the densities of ions 

and electrons are high but substantially equal. It 

embraces the whole space not occupied by sheath. 

The positive columns of the glow and are discharge are ex- 
amples of partially ionized plasmas. Plasma regions are highly 
conducting. In rare gas plasmas, the negative carriers are 
electrons. Ionization of the plasma is produced essentially by 
electron collisions with some photoionization also occurring. 

The positive ions, the electrons and the neutral gas atoms may or 
may not be in thermal equilibrium. The temperature of the positive 
ions is usually greater than the gas temperature, while the elect- 
ron temperature may be very high. Generally, positive ions and 
electrons may be assumed to have a Maxwellian distribution; 
however, variations from this distribution have been noticed 

(oye. 2, 8). 

Even though a plasma, with its ions and electrons moving 
at random according to their characteristic temperatures, is 
essentially neutral when considering a large volume, there must 
be high fields existing at points within the volume. Thus, at 
any point the field will vary with time depending upon the 
instantaneous random configurations of charged particles surround- 
mig Ene point. Fields of this type surrounding a point in a plusma 
may be called microfields and because of them, electrons may gain 


considerable amounts of energy in passing through a series of 


accelerating microfields. 





Another property of an ionized gas which owes its origin to 
Space charge phenomena and high effective conductivity is that of 
electrostatic plasma oscillations. These longitudinal oscillations 
of the plasma are associated with the fact that a small momentary 
displacement of charges of one sign with respect to the other in the 
plasma will establish an electric field (microfield) in such a 
direction as to oppose the displacement. Inertial effects of the 
displaced charges then provide a mechanism for periodic oscillations 
about the point of zero average charge displacement. The frequency 
of these oscillations is determined by the mass of the displaced 
charges and the ‘spring constant' of the electric field produced by 
the displacement (7). 

A simple derivation of longitudinal plasma oscillations might 
be made as follows: 

We write Newton's second law of motion, 

(1) d/dt(m dr/dt) =F =-eE. 

Now considering only wave propagation longitudinally, 
1.e., in the X direction, this reduces to 

(2) mX + eE, =0 

We also write Maxwell's Curl equation 
OE 

gis 


Since no external magnetic field is present, H = 0 


oe Vx Ho= Yt 5 = 


Now reducing the equation to its scalar component in the 
X direction, and denoting partial differentiation with respect 
to time by e , we have 


(4) O=-4 7 neX + Ey, 





Then, integrating and setting the constant of integration 
equal to zero, since we have no external electric field present, 
we obtain 

ae 4 Wt nex 

Substituting (5) into (2), we obtain 

(ye uf no“ X= 0 


This equation gives for the frequency of oscillation 


z, 4 we 
a. f . = ne = 8980 n° cps. 
ti m 
Where: 
Yr is the vector displacement of the electron from 


the point of zero charge displacement 

is the charge of the electron 

is the mass of the electron 

is the restoring force 

rs the eleciric field intensity 

is magnetic field intensity 

is the displacement current density 

is time 

e is the partial derivative with respect to time. 
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Plasma oscillations may be either electronic or ionic in 
nature, and the frequencies are determined by the respective 
electron and ion masses, hence, ion oscillations will be much 
lower in frequency than electron oscillations because of their much 
greater mass (7). 

Plasma electron oscillations may transmit energy by a growing 
wave mechanism. A growing wave results from the excitation of 
plasma electron oscillations and the interaction of these oscillations 
with the primary beam. As the wave moves longitudinally through the 
plasma, it tends to accumulate the slowly drifting groups of oscill- 
ating electrons in advance of the wave and disperse the electron 


groups after passage of the wave. As this redistribution takes place, 


an exchange of energy between electrons and wave in a particular phase 


5 





relationship increases the energy in the wave. The wave amplitude 
grows exponentially reaching the maximum rate of growth when plasma 
frequency approaches the excitation frequency of the wave (8). 

Energy of plasma electron oscillations may also be transmitted 
by fast-moving electrons traversing the oscillating region and having 
their velocities modulated periodically. This type of energy transfer 
obviously depends upon the relative velocities of the oscillating and 
non-oscillating groups of electrons (3). 

Experiments in plasma physics concerned with probe studies and 
plasma oscillations are usually carried out in gas-filled diodes or 
triodes utilizing discharges through rare gases or mercury vapor. 

A diode filled with Argon was used in this study. 

Typically, the voltage profile across such a tube from cathode 
to anode consists of a steep rise in potential near the cathode, 
followed by a slow gradual rise across the positive column of the 
discharge to the anode with a sharp, small rise or fall in the 
region near the anode. The electrons traversing the interelectrode 
distance gain their major potential rise in region adjacent to the 
cathode called the cathode "fall" (9), The cathode is surrounded 
by a sheath of positive ions because of the greater mobility of the 
electrons compared with the ions. Electrons are accelerated as they 
move away from the cathode sheath producing more ionizations and 
attaining a Maxwellian distribution of velocities as they travel 
through the gas making collisions with gas atoms and each other. 

At the negative glow boundary, the entering electrons are composed 
of two or more groups -- fast electrons which Beneuneen produced at 


the cathode and have experienced few collisions in the cathode dark 





space, and slow electrons which have made many inelastic collisions 
in the dark space. The slow electrons have energies bclow the 
ionization potential of the gas, but at, or above, the excitation 
potential; therefore, they make exciting collisions with the gas 
molecules, causing light emission in the negative low region as the 
excited atoms emit quanta in returning to their ground states. 
Afterward, these slow electrons have such low energies that they 
then recombine with positive ions. The group of fast electrons are 
charge carriers which enter into the positive column part of the 
discharge. Some of these electrons may carry through to the anode 
depending upon the mean free path, but more likely, these electrons 
interact in the positive column, give up some of their energy and 
become randomized. The anode will then receive principally elect-~ 
rons coming from the positive column. 

This investigation was originally started at the United States 
Naval Postgraduate School in 1958 by Webb and Garner, continued by 
Flynn and Shoults in 1959, and by Culwell and Johnson in 1960. The 
present work is principally a continuation of the work of Culwell 
and Johnson who redesigned and rebuilt the vacuum system and dis- 
charge tube, but were able to carry out only limited rudimentary 
measurements because of the ee of time. 

Mie principal eLboris of this aoe was in (a) the design and 
development of a Suitable electronic System with which to make 
intensive probe studies in the discharge tube and (b) then the 
making of such measurements. A unique system was designed yielding 


Continuous graphical plots of conventional Langmuir data such as 





probe current vs. applied probe potential and probe current 
squared vs. applied probe potential. This system shows promise 
of extension to yield accurate space potential determination and 
electron energy distribution functions. 

This report covers plasma oscillation and probe studies made 


in low pressure argon discharges. 





2. Previous Work 

Investigation of gaseous discharge phenomena at low pressures 
was given its first real impetus in 1924 when Langmuir and Mott- 
Smith (10) reported a new method of making extensive studies of 
electrical discharges in gases at low pressures. Their technique 
consisted of placing a small auxiliary electrode (commonly referred 
to as a probe) in the discharge” and studying the complete probe 
volt-ampere characteristic using the collector theory proposed 
by Langmuir (11). 

Further studies by Langmuir (12) in 1925 disclosed that appar- 
ently there were three energy groups of electrons present in the 
discharge. He classified these groups as "primary electrons" 
(those electrons which keep nearly all the momentum gained in the 
cathode ion sheath), “secondary electrons” (those electrons having 
a Maxwellian distribution of energies proportional to the energy of 
the primaries), and a third group of “ultimate electrons" (those 
electrons having a Maxwellian distribution of energies independent 
of the primaries and much lower than the secondaries) . 

These observations indicated that there must be some type of 
energy transfer mechanism involved to permit certain electrons to 
gain energies in excess of the tube drop while others became energy 


deficient. Langmuir referred to this phenomenon as scattering. 


"Floating probes had been used previously in analyzing discharges, 
however, they were of limited value because of the failure to appre- 
Ciate sheath formation surrounding the probe and the consequent 
failure to measure space potential. 





Dittmer (13) experimentally disproved the hypothesis that a 
new mechanism (such as the possible large increase in the target 
area of excited molecules) was the cause of scattering by showing 
that the scattering reached a constant value at some distance from 
the cathode. Dittmer attempted to find oscillations in an effort 
to validate the hypothesis that electrical oscillations in the dis- 
charge caused scattering. He was unsuccessful in all attempts at 
minding such oscillations. 

Tonks and Langmuir (2) developed a simple theory for electron 
and ion oscillations which was analogous to the electrons being im- 
bedded in a jelly-like substance composed of ions. The electrons 


have a natural frequency given by 


ne = ne“ cps 
TT mM, 


where n is the electron density per cubic centimeter, e the electronic 


nv" 





eeerse. mo the mass of the electron. The ions obey a similar law with 


a natural frequency approaching 


BF aie sl (me/m,)* 
where ms is the mass of the ion. 

Tentative verification of this theory was obtained by detecting 
oscillations in the 300 to 1000 megacycle range which corresponded 
reasonably well with the predicted frequencies. It was conjectured 
that additional oscillations in the 50 to 60 megacycle, 1.5 megacycle 


and lower ranges might correspond to oscillations of the primary 


beam electrons and ions, respectively. 
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Merrill and Webb (14) carried on extensive investigations in 
a mercury arc discharge at low pressures and discharge currents 
with a movable probe to determine a possible relation between the 
plasma oscillations and scattering as suggested by Tonks and 
Langmuir. They found that scattering occurred in several narrow 
regions and was associated with periodic oscillations of consid- 
erable magnitude which agreed with the Langmuir frequency within 
ten per cent. The fact that several regions were observed has 
been explained by Bohm and Gross to be a bunching effect of the 
electron beam. Merrill and Webb further noted that for each 
scattering region there was only one observable frequency and 
that the maximum intensity of the oscillations occurred at 
from one half to one millimeter from the scattering region on 
the anode side. 

Merrill and Webb concluded that the scattering was due to 
the oscillations receiving energy from the fast electrons, which 
then give some of the electrons excess energies, and that the 
observed difference in the positions of the oscillations and the 


* 
scattering was a consequence of the experimental technique. 


“With the recent theory proposed by Allis and partially 
confirmed by Emeleus, one might find that an analysis of Merrill 
and Webb's 600 complete volt-ampere characteristic curves would 
show that their observation between scattering and oscillation 
locations was of considerable physical significance. 


ileal 





Although no major breakthrough in understanding oscillations 
occurred in the next ten years, many important general observations 
were made. A few of these are stated below (15). 

(a) The operating position on the current - voltage char- 
acteristic of the tube under investigation may affect 
the amplitude of oscillations observable in external 
circuits. It has been found that oscillations can 
usually be detected over wide sections of the 
characteristic for oxide-coated cathodes. 

(b) The position of the probe in the discharge markedly 
affects the amplitude of the observed oscillations; 
however, the observed frequency usually remains 
about the same. 

(c) The signal amplitude of the oscillations depends 
considerably on the electrode geometry and upon 
which electrodes are used for detection of 
oscillations. 

Emeleus (16) has pointed out the value of evel observations 
in electrical discharge studies using oxide-coated cathodes in ar- 
gon at a few microns of pressure with a low voltage discharge. A 
representative sketch of a gaseous discharge in which high frequen- 


cy oscillations are present may be shown as follows: 


: 1 2 3 y A 
/ — 
5 
5 
i. —_—a. 


Fig. A. C, cathode; 1, cathode sheath; 2, unscattered beam 
section of discharge; 3, meniscus; 4,5, scattered beams; A, anode. 
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The principal static fall in potential occurs in region lL; 
the rest of the discharge is assentially at anode potential. 
Between 1 and 2 represents the familiar Dittmer scattering region 
with dispersion occurring in or just after 3. It has further been 
noted that strong oscillations have only been detected to the 
Moat OL 3. 

Prior to 1950 little attention had been given to the effect 
of probe contamination on the validity of the volt-ampere probe 
characteristic. Easley (17) in probe studies of a mercury dis- 
charge buffered with argon showed, by drawing sufficient electron 
current to make the probe cherry-red (and thus outgas it by ohmic 
heating) before each reading that many of the non-linear logarith- 
mic plots previously obtained were due to a change in contact 
potential and thus did not indicate the presence of a non-Maxwell- 
ian distribution of electron velocities, but rather probe contam- 
ination. 

The results of Oleson and Found (18) have shown that non- 
linear logarithmic plots were obtained in discharges in which 
moving striations were present. The logarithmic plots obtained 
looked as though two Maxwellian distributions were present in 
the discharge. However, it had become fairly evident that in 
many cases (such as in discharges at pressures of the order of 
several microns of mercury) that a Maxwellian distribution of 
energies did not exist. Consequently, increased effort was 
directed toward new techniques of determining the actual energy 


distribution function rather than assuming it to be Maxwellian. 


13 












===> * 
Guam = ee » Ot = ih 
—_——_« «6. 
- a 1 
* —_ 
pau gga 
7 Sa es o—- <—-_> 1 





-_— —x —— 
a a 


<i fl 









a 


—_ 


eaa— 


It is worthwhile to note that Langmuir considered the energy 
distribution function problem in his original theory on collectors, 
and though he did not employ the technique, he did show that ob- 
taining the second derivative of probe current with respect to 
probe voltage would give the true energy distuinatene: 

Haigh and Smith (19) undertook a study of the energy distri- 
butions in the positive columns of low pressure (.005 mm Hg + 
.085 mm Hg), hot cathode discharges in argon and air by a magnetic 
spectrum method. A continuation of this study by Haigh and 
McCulloch (20) obtained data on air, argon and nitrogen. Their 


results, using argon, are shown in Figs. B through E. These 


plots illustrate the types of distributions obtained in argon. 
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“Strictly speaking, this is true for a cylindrical collector if 
the velocity distribution is isotropic. In cases in which one 
doubts the presence of an isotropic distribution, a spherical 
collector is required. This is shown in Langmuir's original 
development. 
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it is interesting to note that for a given pressure, a decrease 
in current has the same effect as a decrease in pressure at constant 
current. Haigh and McCulloch hypothesized that the second peak in 
their distributions may be the result of the "combined effects of the 
Ramsauer free paths and the tendency for a larger proportion of high 
velocity electrons as the pressure decreases." 


As far back as 1939 Druyvesteyn (21) had shown that the electron 


energy distribution fun¢ction could be found from the expression 


- 
om V : a°t, 
N f(energy) =N#(eV) = 2 Se a 
‘Se ea Sea 
dV, 


where N@(eV) is the number of particles per unit of energy per unit 
volume, and in this work, customarily put in units of reciprocal 


electron volts and reciprocal cubic centimeters, i.e., ev! om-?: 
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A is the area of the probe, Mm. is the mass of the clectron, 


is the probe current and Vo 1 


the probe voltage. All are expressed in consistent units. 


e is the eleetreamic charge, iy 

Despite this, no significant effort was expended to obtain 
energy distributions because of the inherent inaccuracies in- 
volved in the tedious graphical method of obtaining the second 
derivatives of probe current with respect to probe voltage. 

Sloane and MacGregor (22) developed an idea originally pro- 
posed by Emeleus for electrically obtaining the energy distribution. 
By putting a small alternating signal on the probe it was shown 
that the de increase caused by the ac signal was directly propor- 
tional to the second derivative of probe current with respect to 
probe voltage. Though it appeared that this technique might yield 
the first real gain in obtaining distribution functions, Emeleus 
mentioned in Loeb (3) that "not too much confidence should be 
placed in the distribution derived by this technique ey 
principally because of the many ways the plasma might react to the 
alternating potential. 

Recently, Boyd and Twiddy (6) used an improved version of the 
Sloane and Mac Gregor technique by amplitude modulating that ac 


signal thus permitting work in plasma having higher noise levels. 
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A typical energy distribution measurement by Boyd and Twiddy 
for a striated head in hydrogen shows at least two electron energy 


Bremps, aS indicated in Pig. F. 
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Electron Energy (eV) 
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Boyd and Twiddy have also shown that for a non-isotropic, 
strongly directed discharge a cylindrical probe will give an energy 
distribution which will have negative values of energy between two 
major electron groups. In every case studied by Boyd and Twiddy, 
two well-separated groups of electrons (and sometimes a third group) 
were detected. Further, it has been suggested that the energy ex- 
change between the electron groups might be a much more rare phenom- 


enon than had been expected. ie 





In 1956, Medicus (23) developed a simple technique for deter- 
mining the energy distribution function from the probe character- 
istic curves plotted automatically with an X - Y recorder. 
Medicus' method is a graphical one in which tangents are drawn to 
the probe curves and the differences between the tangents repre- 
sent the second derivatives. Fig. G below outlines the method 


pieterially, 


1=0 
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To obtain the best results the probe curve should be ampli- 
fied so that all regions of interest lie at about a 45° slope. 
The actual probe curves may be used directly only if the ion current 
is small compared to the electron current, or if the ion current 
varies linearly with voltage; otherwise, one must use the net 


electron current. 
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Medicus (24) obtained results which seem to further substanti- 
ate the suggestion of Boyd and Twiddy concerning the possible low 
energy exchange between slow and fast energy groups. Medicus showed 
that when a discharge consists of both a high energy, non-Maxwellian 
group of electrons and a low energy Maxwellian group that "the 
characteristic linear part of the corresponding probe curves, extra- 


polated to plasma potential, directly determines the random current 





density of the fast group." The expression for this relationship 
is given by 4 1p 
Sp 
eA 


Fig. H shows how this is obtained graphically and how it is re- 


lated to the high energy electron distribution. 
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Applying the diffusion law to his result, Medicus found that 
the interaction between the slow and fast groups must be small and 
that cases may exist “at low total plasma densities where the 
particle density of the slow group is always comparable to, or 
even much higher than, that of the fast group." . 

Continuing the search to improve probe techniques considerable 
effort has been expended to develop efficient methods of eliminating 
probe contamination while obtaining the probe characteristic curves. 

Wehner and Medicus (25), using an automatic X - Y plotter, 
were able to obtain a complete probe characteristic Curve in sever- 
al seconds coming from a cherry-red probe condition. They claimed 
no deviation from a Maxwellian distribution over nearly 4 powers 
of 10 in probe current even in tubes with oxide-coated cathodes, 
providing the characteristic was recorded in several seconds after 
the probe had been thoroughly cleaned as described above. 

White (26), while making measurements in a hollow-cathode 
glow discharge, found that excessive material was sputtering from 
the cathode onto the probe. It was found, that by keeping the 
probe at 30 volts below floating potential and by pulsing it 
every 0.1 seconds for 0.004 seconds, it was possible to maintain 
a clean probe by ion bombardment and hence, obtain a reliable 


probe characteristic. 


* : ; : a ae : 
This would imply that the particle densities of the different 
energy groups would have discrete values at any point in the 
discharge tube for a given set of parameters. 
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several months later Waymouth (27) reported a Ssimiiar 


pulse technique for obtaining "uncontaminated" probe 


characteristics. 


Some results of electron current vs. probe voltage for 


discharges at pressures of six microns of mercury plus one 


mn of neon are shown in Fig. I, below, taken from Waymouth's 


avercle (279% 
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An attempt to obtain energy distribution functions by the 
method of Medicus gave poor results because of the small size of 


the oscilloscope photographs. Standard scatter of data is shawn 


for a probe 0.25 cm from the cathode. 
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It was found that this pulse technique” represents at 
least a 600-fold reduction in the net rate of probe contum-~ 
ination in gas discharges. It was also shown that the energy 
distribution of electrons in the negative glow consists of 
two groups, a Maxwellian group and a high energy group with a 


potential approximately equal to the cathode fall. 


“Ie Re important that in the zeal to overcome one problem 
it is not replaced with a more complex one. For these measure- 
ments to have meaning the ion sheath must be in equilibrium; 
therefore, it is clear that there is some upper limit on the 
rise time of the pulse so that the sheath thickness is essenti- 
ally in constant equilibrium, i.e., where the distance the ions 
Can travel during a pulse equals the product of the pulse dura- 
tion, the ion mobility and the sheath field. 
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Waymouth (28) has pointed out an important new 
diagnostic tool in making a rapid evaluation of tne 
principal electron energy groups directly from the Linear 
iy vs Vp probe characteristic curves. His technique is 


essentially this. 


ee & 
3 < 
N@ (eV) ~ Vp d i,p/dv, : 


therefore, whenever the linear 1p vs Vy curve has a constant 
slope 
N@ (eV) = 0, 
i.e., there are no elcetrens On corresponding energy in 
the energy distribution spectrum. 


An example is shown in the following plots: 
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A modern radio frequency technique a determining the 
plasma resonant frequency in a gaseous discharge has been 
developed by Takayama, Ihegumi, and Miyozaki (29). They 
discovered that if an rf signal is superposed on a probe, 
which is maintained negative at all times rclative to space 
potential, then a resonant current increase appears in the 
circuit at the Langmuir frequency, f = (pes )" 


cps 
TT Me 
The relation between the resonance and the frequency of 


imposed rf voltage is as shown in Fig. L, below, for a 


homogeneous mercury plasma at Q.1 mm Hg pressure. 
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Fig. bL. 


Takayama, et al., suggest that this is a beam-plasma 
type of interaction and that the resonance effect will not 
be discernible when the energy supplied by the rf field is 
comparable to the thermal energy of the electrons. They 
point out that an increase in pressure will likewise causc a 
decrease in resonance since the height of the resonance peak 
also depends upon the collision frequency of the electrons with 


atoms. ou 





In 1957, Allis (30) suggested that as an electron beam 


travels through a region of increasing plasma density, electron 


oscillations are being generated in 
not propagate in the absence of the 
an electron beam travelling through 
density might generate oscillations 
outrun the beam. 

Allis further pointed out that 
electrons is going through a plasma 
is reduced by the growing amplitude 


energy being transferred into field 


a region in which they can- 
beam. On the other hand, 
a plasma of decreasing 
which could propagate and 

. 
if a monoenergetic beam of 
then the average beam velocity 
of the oscillations, the 


energy in the wave and in- 


creased kinetic energy of the random plasma electron group. 


Allis states, "growth occurs mainly 


below the plasma frequency" 


So that the free plasma cannot transmit these waves, {see Fig. M). 


He further mentions that if the beam should peter out through 


scattering or slowing down, the plasma waves must decay expon- 


entially along with the beam. 
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Pig. M. Amplification factor for an electron 
beam in a plasma (30). 
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A recent investigation of electron beam interaction with a 
plasma was conducted by a group of Russian workers (31). In their 
paper they point out that theoretical investigations have shown 
that the state of an unmodulated beam of electrons moving through 
a homogeneous plasma is unstable. Density and velocity fluctuations 
arising in the beam or the plasma, may excite longtudinal electro- 
magnetic waves (plasma oscillations) of growing amplitude in the 
plasma. They note that the initial fluctuations may be caused by 
external forces. 

These workers state that the observed maximum excited or 


amplified frequency is determined by the following relation: 


be 


Wo = 2p ony S 





nt 





Where, 

“5 is the observed frequency in radians 

Wp is the Lnagmuir plasma frequency 

Vis the beam velocity 

S is the average thermal plasma velocity 
They also point out that a large velocity spread in the beam 

may remove the instability. If the regularly moving beam consists 
Gf ions, the excitatiom is less, not by a factor Tia7 ly eee Un we nidny 
aa Therefore, the interaction mechani#sm considered is 


also important if an ion beam moves through a plasma. 
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Plasma instabilities may be either convective or absolute. 
Convective instabilities tend to move with the flow and the 
oscillations they produce amplify or diminish with time. Ab- 
solute instabilities produce oscillations which grow contin- 
uouSly with time. 

Emeleus and Garscadden (32) have very recently done some 
work in exploring plasma oscillations with a transverse electron 
beam observing the resultant energy spread of the transverse beam 
when if iS in, or has traversed, a section of the main discharge 
where scattering or oscillations can be detected. 

A typical result, for a discharge of 1.5 microns pressure 
in mercury vapor with a discharge current of 27 ma (8 mm diameter 
cathode) with an electrode separation of 42 mm and a meniscus 


forming LO mm from the cathode, is shown below. 
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Fig. N. Energy spread for a transverse electron 
beam in a meniscus-type discharge. 
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The conclusion is made that “with the discharge described 
* * * The results appear to be compatible with considerable 
randomization and possible break-up of the main beam by the 
oscillation" as the beam interacts with the plasma in the 

. . * 
meniscus region. 

Very recent investigations by Mahaffey, McCullagh, Garscadden 
and Emeleus (33) are pertinent in correlating growth or decay of 
oscillations with plasma density gradients. 

Their work was also done in a mercury discharge at approxi- 
mately one micron of pressure using an eight mm oxide-coated 
cathode with a tube current of 27 ma and a tube voltage of about 
20 volts. Two main sets of measurements were made with emphasis 
on (a) frequencies, intensities, and energy spread of the electron 
beam (found from the break in the ip vs. Vy) curve at which fast 
electrons began to be received). (b) frequency measurements, 
electron temperatures, amd electron concentrations (found from 


a 


p Vs Vp curves near space potential). 


They conclude that: 


Some correlation has been found between growth or 
decay of oscillations and motion of the electron 
beam down or up a plasma concentration gradient, 
which appears to be in general accord with the 
plasma dispersion theory. 


*of. our results, Section 5 of this report, for a beam generated 
in the discharge tube interacting with the plasma. 
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3. Equipment. 


The equipment used can be divided into four main component 
groupings as follows: 

(1) Gaseous discharge tube in which a glow discharge could 
be continuously maintained. 

(2) Vacuum system for evacuation and filling of the tube. 

(3) Frequency detection and measuring equipment. 

(4) de electrical circuit for supplying power and for 
controlling and measuring discharge characteristics. 

Figs. 1 through 3 are photographs of the general experi- 
mental setup. 

a. Discharge Tube. Figs. Ya and 4 are photographs and Fig. 
5 is a detailed drawing of the tube (diode). This is the same 
tube used by Culwell and Johnson in their work during 1960. No 
changes were made to the tube. The electrodes and probes were 
attached to soft iron cylindrical cores sealed in glass which 
made possible the variation of interelectrode distances and 
probe positions by moving the iron cores with externally 
applied magnetic fields. 

The probe used in this work was a cylindrical probe made 
of .008 inch diameter tungsten wire enclosed in a glass sheath 
Mpeteease millimeters from the tip. 

The ring probe shown in Fig. 5 was not used during the 
experiment. It was damaged during electrode movement, breaking 
loose from its mounting and falling to the bottom of the discharge 


tube. 
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The cathode was a Philips impregnated "L" cathode with a 
two centimeter diameter face indirectly heated by a tungsten 
element drawing seven amperes dc. 

The anode was a three centimeter diameter nickel disc. 

This tube was originally designed with high purity of the 
gas being a primary design criterion. As such, the coiled wire 
leads from the electrodes and the probe to the external connect- 
ions were unavoidable. The inductive effects of these leads may 
be considerable at electron oscillation frequencies and tend to 
mitigate against making large amplitude high frequency measure- 
ments through them. 

b. Vacuum System. The vacuum system built by Culwell and 
Johnson (8) was not changed. It remained operative during the 
entire course of this work without leaking. Pressures of the 


2 mm of mercury were readily attained. All valves 


order of 10° 
in the system were glass stopcocks lubricated with Type N Apiezon 
grease. The diffusion pump was a three-stage air-cooled pump 
using Octoil-S. The major part of the system was baked out with 
a portable oven (see Fig. 1 ). Heating tapes were used to bake 
out external glasstubing. An induction heater was used to out- 


gas the anode. Oyutgassing of the Langmuir probe was accomplished 


by ohmic heating while the tube was operating. 


“Manufactured by Semicon of California, Inc. 
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Low pressures were meaSured with a Veeco Type RG 31 A 
ionization gage. Pressures in the operating range of one to 
ten microns were measured with a Langmuir viscosity gage (35). 
Tube voltage readings were also used as an indication of tube 
pressures after the voltage-pressure relation had been determined. 
Tube voltage is a function only of tube pressure if other para- 
meters such as anode-cathode spacing and discharge current remain 
constant; therefore, a given pressure could be readily duplicated 
by evacuating or filling the tube until the correct voltage was 
obtained. Thus, when fresh argon was desired in the tube, flush- 
ing and refilling was quickly and accurately accomplished without 
relying on the Slower pressure measurement process using the 
viscosity gage. 

c. Frequency Detection Equipment. A cavity wavemeter, Test 
Set, Model AN/UPM - 2, with a frequency range of 80 to 1220 Mcs 
was used to detect electron oscillation frequencies by coupling 
with the wire probe at anode potential. The dc signal from the 
cavity wavemeter detector crystal was fed into a Hewlett-Packard 
Model 425 A de micro volt-ammeter to obtain detectable signals. 
Currents as small as 1071¢ amperes could be detected with this 
instrument. 

A Hewlett-Packard UHF signal generator, Model 612 A, with 
a frequency range of 450 to 1230 Mcs was also tried utilizing the 
technique of parametric amplification suggested by K. Takayama, 


et. al. (29) which essentially consisted of putting the rf 


ue 





potential on the probe held at space potential while the frequency 
of the generator was varied in an attempt to cause resonance at 
the electron oscillation frequency. 

Ion oscillations were detected with a Tektronix Type 581 
oscilloscope and also simultaneously counted with a Hewlett-Packard 
Electronic Counter, Model 521 A. The ion signals were amplified 
by a Hewlett-Packard 40 db amplifier, Model 450 A, to feed count- 
able signals to the counter. 

d. de Circuitry Components. Fig. 6 is a schematic block dia- 
gram showing the general de Circuitry employed. 

Filament voltage was supplied by storage batteries regulated .. 
by a charging circuit during operation. Seven amperes dc was 
Sufficient to heat the Philips L type cathode to 1050° C. Cathode 
temperature was measured on the face with a Leeds and Northrup 
optical pyrometer. Cathode temperature was maintained within + 
SC during operation. 

Electrode spacing was measured with a locally manufactured 
cathetometer which allowed electrode positioning within + 0.25 mm. 

The discharge tube was powered by a Kepco Labs Regulated 
High Voltage Power Supply, Model 1250 B, having a voltage range 
from 0 to 1000 volts and a current output from O to 500 milliam- 
peres. 

An F. W. Moseley Autograf X - Y Recorder, Model 2, was used 
in the plate circuit to obtain continuous plots of the volt - 
ampere tube characteristic simultaneously with the taking of probe 
data, thus the effect on discharge current of drawing probe 


current could be monitored and recorded. 
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tui@ial at®eMts at takiine probe data were designed around 
osealogwa@phicerecording of the probe characteristic curves. 
Spurious high frequency signals in the range 0.5 to 5 Mcs of un- 
known origin interfered with these attempts despite coaxial 
Shielding, etc. The oscillographic method was therefore abandoned 
aS a primary means of meaSurement and a system designed to mini- 
mize ac interference was constructed. A description of this system 
follows. 

The probe voltage was supplied by a Hewlett - Packard Low 
Frequency Function Generator, Model 202 A, set for triangle wave 
generation to obtain a voltage source varying linearly with time. 
This signal was then amplified and biased to obtain the desired 
positive and negative voltage excursions relative to the anode. 

A Donner Analog Computer, Model 30, with Problem Board, Model 
30 - 6, was used to amplify the voltage drop across the variable 
Besistor in thesprobe circuit, thus, producing a voltage proportion- 
al to probe current in order to ensure a voltage range suitable for 
best operation of the function multiplier and X - Y recorder. Ten 
operational amplifiers were available with this computer. The 
analog computer was also used as a source of five separate regulat- 
ed dec power supplies from its initial conditions panel. The analog 
computer amplifiers operate in the range -100 to +100 volts with 
long time stability better than 0.5% full scale value. Initial 
condition voltage sources operate from 0 to 100 volts. Computer 
amplifiers have a maximum drift characteristic of 4 millivo)ts 


fer hour, with a short term random drift less than + 2 mv. 
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Balancing of amplifiers is easily accomplished and holds drift 
error” temneeeeole levels over short time intervals for high 
Signal to drift ratios. Potentiometers on the problem board 
allowed precise gain control on amplified signals. 

A Donner Electronic Function Multiplier, Model 33, was used 
to square the values of voltage proportional to probe current. 
iis eQuipiem: would accept inputs of + 100 volts in either or both 
input channels (X, Y) and produced an output voltage of -0.01 XY. 

An F. W. Moseley Autograf X - Y Recorder, Model 2S, was used 
to make continuous plots of probe data. Probe voltage with re- 
spect to the anode was applied to the horizontal axis which has 
a range from O-7 .5 millivolts to O-150 volts full scale. The 
voltages proportional to probe current, or probe current squared, 
were applied to the vertical axis which has a range 0-5 millivolts 
to 0 - 100 volts full scale. Input impedances of both axes are of 
the order of 2 megohms. Sweep rate was controlled by the low 
frequency function generator. It was generally kept from 0.05 
eos to O.2 eps to @v0id inertial error in the X - Y recorder. 
This error is negligible at these sweep rates. The Autograf is 
accurate to 0.25% of full scale value on all ranges when calibrated. 
PaLibration dritt is negligible. 

Both the output of the function multiplier and input of the 
X - Y recorder have low pass filters keeping ac interference to 


ani ni mim . 
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A Lambda Electronic Corporation Regulated Power Supply, 
Model 71, was used as a probe voltage supply when a non time- 
varying source was desired. 

A Donner Variable Base Function Generator, Model 3750, was 
also tried in an attempt to automatically make plots of the 
natural logarithm of probe current vs. applied probe potential. 
Preliminary attempts were unsuccessful because of difficulty in 
programming a suitable logarithmic function. This technique was 
temporarily abandoned because of the press of time. It does 
appear quite feasible, nevertheless, and pursuit along this line 


* 
might prove to be quite rewarding. 


+ 
F. W. Moseley Co., manufactures a logarithmic attachment for 

the X - Y recorder which enables semi-logarithmic plots to be 

made without need of a function generator. They also manufacture 

a curve follower attachment to the X - Y recorder which allows 

functions to be drawn on paper with magnetic ink; the pick-up 

coil will then feed out the function desired as the independent 

voltage is driven onto the horizontal axis. A system such as ours, 


utilizing this technique, would require two recorders -- one to 
record, one to generate functions -- both being driven in time 
synchronism. 
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WU. Procedure. 

Prior to making a data run, the vacuum system was evacuated 
to as low a level as possible -- in the range of ae mn of mercury. 
Then the system was baked out with the oven at 450° C for a period 
of between twenty four and thirty six hours. During this time the 
tubing and liquid air traps external to the oven was wrapped in 
heating tape and heated to 350° C. Following bake out, the cathode 
was heated to operating temperature with the filament and the anode 
was Outgassed with the induction Hosa The system was then 


-8 and 107° 


allowed to pump down to between 5 x 10 mm of mercury. 
This normally took 24 hours after completion of bake out. The 
system was then filled with approximately 2 mm of spectroscopic- 
ally pure argon. This filling procedure was repeated several. 
times in order to flush the tube with the inert gas. With argon 
in the system, the tube could be evacuated down to 5 x 1072 mm 
of mercury in about one and one half hours. The final pressure 
was obtained by discharging the tube, measuring the pressure with 
the Langmuir viscosity gage, and correlating it with the plate 
voltage; then, evacuating the tube in a small increment of pressure, 
while the tube was still discharging, the process of measuring 
pressure and correlating with tube voltage was repeated until 
pressure ~ tube voltage calibration had been achieved. The final 
pressure was then readily obtained by evacuating until the proper 
plate voltage was obtained. 

“This bake out and purification procedure is tedious, but necess- 
ary. For example, failure to be careful about extreme cleanliness 
results in a pulsating discharge with a frequency of one to 25 cps 


and an extremely unstable tube characteristic (i.e., the tube "runs 
away" at very low discharge currents). 
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An arbitrary interelectrode separation distance of about 42 
mn was selected so that probe measurements could be made with the 
probe at increments of two mm along the axial length of the tube 
between electrodes. Thus, about twenty different positions between 
anode and cathode were studied at approximately equal intervals. 

Electrode positions were calculated and set at the catheto- 
meter as follows: The cathetometer was set at its reference value 
for the cathode. The cathode was then moved using the magnetic 
coil over the glass surrounding the soft iron magnet within the 
tube until the cathode was lined up with the reference mark on the 
cathetometer. The cathetometer was then traversed to the probe 
position where its value was recorded (since the probe position 
remained fixed, this position was used to estimate the statistical 
accuracy of the cathetometer measuring system). The cathetometer 
was then traversed to the anode position and the anode moved to 
coincide with the reference mark as described for the cathode. 

All cathetometer peeimes were made going in the same direction 
to minimize backlash errors. 

A tube volt-ampere characteristic trace was then made and 
recorded on the plate circuit X - Y recorder to determine the 
operating range of tube currents. When the gas in the tube was 
very clean the tube remained stable up to 200 milliamperes dis- 
charge current. Over this range, then, a group of representative 
discharge currents were selected to study which gave suitable 
values of probe current (i.e., probe currents which did not pro- 
duce voltages in excess of 100 volts to the operational amplifiers 


of the analog computer). The discharge currents selected were 


Si 





25, 50, 75, 100 and 125 milliamperes. All of these currents were 
attainable using the voltage control on the Kepceo Power Supply, 
the plate circuit resistance remaining fixed. 

Proceeding next to the probe circuit, calibration procedures 
were employed to minimize all electronic errors. The amplifiers 
of the analog computer were balanced to minimize drift. The 
function multiplier was balanced for drift and accurate suppression 
of cross products. A potentiometer on the problem board was used 
to calibrate the maximum peak-to-peak excursion of the triangle 
probe voltage supply to an absolute magnitude of 150 volts as 
traced on the horizontal axis of the X - Y recorder. This voltage 
was then biased with an initial condition voltage source (set with 
a Donner Null VTIVM) to give the maximum positive and negative ex- 
cursions desired. Those used were +20 volts to -130 volts with re- 
spect to the anode. 

Calibration of the vertical axis voltage corresponding to the 
probe current was accomplished with a problem board potentiometer 
after the signal had been amplified in the two nominal gain-ten 
stages. This potentiometer was then adjusted to give the same 
magnitude at penpoint on the X - Y recorder as the ammeter reading 
probe current. This had to be accomplished dynamically, but could 
be done to quite good precision (as well as one could read the 
meter) at very slow sweep speed. This was commonly done at the 
maXimum ammeter and pen excursion when both momentarily came to 


rest. Inertial effects were considered to be negligible. 
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Once these calibration procedures had been made, data runs 
were started. For each position of the probe a linear probe 
characteristic curve was made for each of the five tube discharge 
currents cited previously. Blowup traces of these curves were 
also made of the ion (negative) portions of the trace. The paper 
was then changed in the X - Y recorder and the current squared vs 
applied probe potential plots were made for each of the discharge 
currents. Diiferent vertical scales were used for electron and 
ion regions of the curves corresponding to those used for the a 
linear traces. Simultaneous with these measurements, the plate 
circuit volt-ampere characteristic was being monitored and traced 
giving a record of its changes during the probe characteristic run. 

Prior to each new position of the probe the calibration 
procedures cited above were employed, except that analog computer 
amplifiers were only balanced approximately hourly. 

A run normally took twenty to twenty five minutes to make -- 
about five minutes to trace the curves, and about fifteen to twenty 
minutes to reposition electrodes and calibrate equipment. Twenty 
such positions between anode and cathode were studied using this 
procedure during a twelve hour period yielding four hundred 
separate probe characteristic traces. 

Reproducibility of curves is quite good within the limit of 


electronic errors which will be commented upon later in this report. 
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A subsequent, more limited, data run was made for ten equi- 
spaced probe positions (electrode separation 40 mm, argon pressure 
6.9 microns of mercury, discharge current 50 ma), for the purpose 
of obtaining second derivatives of probe current with respect to 
probe voltage. Two simple RC differentiating circuits in cascade 
were used to do this on the probe current signal, (see Fig. 6). 

Thse derivative curves were used then to obtain electron 
energy distribution functions. 

On this data run, ion oscillations were measured by 
oscilloscopic photographs and electronic counter measurements. 
Electron oscillations were measured with the cavity wavemeter 
coupled to the dc micro-ammeter. 

During this run great care was taken to evaluate the effect 
of probe contamination by retracing a curve after various time 
intervals and noting trace deviations from the original; then 
after cleaning the probe by ohmic heating, the curve was again 
retraced to note the tendency to return toward the original 
trace. The results of these tests were quite conclusive and 


are stated later in Section 5, Results. 
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5. Results. 

Plots of two representative probe positions from the first 
run (consisting of twenty probe positions and five discharge 
currents) are shown in Figs. 7 through 13 for the following 
func tional relationships: 


2 


iy me, V5 1 vs V , Ini. vs V_, Ig vs Vg, where the sub- 
oer Pp p 


p 
script p denotes the probe parameter, subscript + denotes the ion 
region of the probe characteristic curve, subscript - denotes the 
electron region of the probe characteristic curve, and subscript 
d denotes the tube discharge (or plate) parameter. 

A summary of three important plasma parameters for all 
twenty probe positions between anode and cathode are shown in 
plots as indicated below: 

(a) n_ (electron density) vs X (cathode - probe distance) 
obtained by using the slope of bg vs V5 curves (Langmuir Theory) 
is shown in Fig. LH. 

(b) Vsp (space potential) vs X, obtained from ln ip vs Vp 
plots at the first positive voltage deviation from linearity is 
SvoOwn im Fig. 15. 

(c) n.vs J (discharge current density) for ten probe 
positions spaced approximately 4 mm apart to correspond with the 
positions used during the second data run is shown in Fig. 16. 
These data are perhaps capable of being correlated with the re- 
sults of the second data run, which was made for only one discharge 


current (50 ma), despite the fact that the two runs were made at 


slightly different pressures. 
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A second data run was made at 50 ma discharge current for 
10 probe positions spaced at4¥mm intervals with an anode - cathode 
separation distance of 40 mm. The gas pressure was 6.9 microns of 
Hg. Results of this run are indicated below: 

(a) Plots of ip vs Vp, di,/qv,, : d°ip/ave are shown in 
Figs. 17 a, b through 18 a, b for two representative probe posi- 
tions. 

(b) Electron energy distribution functions for all ten probe 
positions are shown in Figs. 19 through 28. These functions were 


obtained from the second derivative of probe current with respect 


to probe voltage according to the relation 


eC peer oe 


e 
avi, 


NB(eV) = 2 20%» | pn 


Where N@(eV) is the number of electrons per electron volt of energy 
per cubic centimeter, A is the area of the probe, m is the mass of 
the electron, e is the electronic charge. For the geometry of the 


probe used in this work, the function reduces to 


nn 


iL3 2 





NB(eV) = 1.722 x10 Vi, ai, tar 
ave 
p 


Two things are to be noted concerning these derivatives and 
associated energy distribution functions. First, the differentiat- 
ing circuit was not perfect -- in particular it appears that the 
Capacitor was not unloading quite fast enough to give tne values 
when the derivative was decreasing. This error was relatively 


small and was arbitrarily adjusted to make it fit the first 
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derivative. Second, the X - Y plotter developed some mechanical 
troubles during this run and gave a magnitude for the second 
derivative of approximately 1/3 of the true value. Approximate 
calibration of the derivatives was achieved by using another 
machine which was reading true and also numerically differentiating 
(well realizing the errors inherent in this process) the in VS ‘ 
curve to establish the scale factor. The second derivatives and 
the distribution functions shown in Figs. 19 through 28 are shown 


in their original scale as they came from the machine. The 


correction factor is applied to the constant term, i.e. K 





correction 
13 di 
coer ee x TO Me Pp . This correction factor may still be 
ave 
Pp 


in error by perhaps ten percent; however, these results are shown 
because the technique appears valid and the information the curves 
yield appears to be a valuable diagnostic tool. The actual magni- 
tudes really are not critical for diagnostic purposes; shapes and 
changes in shape in going from one probe location to another are 
important because these factors show how the energy of the particles 
is being transferred. 

(c) A plot of n_ vs X is shown in Fig. 29. These values 
were obtained by integrating the energy distribution functions. By 
this method the density could be divided between the high and low 
energy groups and changes in density groupings noted as the probe 


position was changed. The integration was carried out numerically 
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using the trapezoidal rule.” Results were checked with a plani- 
meter and agree within three per cent. Fig. 4/7 shows a typical 
distribution function, and the technique employed to integrate it 
and obtain the various densities. 

(dq) A plot of n_ vs X is shown in Fig. 30. These data were 


obtained by the conventional Langmuir method of using the slope of 


2 
i 

Pp: 
plots have the same general shapes, and essentially the same peaks 


Vs Vy plots. Comparison of Fig. 29 and Fig. 30 shows that both 


and gradients are obtained by both methods. The absolute magni- 
tudes vary by as much as 30%. This may be attributable to the 

scale factor assumed for the second derivative method; but, possibly 
part of the difference (perhaps ten per cent) may exist in comput- 
ing the slopes by the Langmuir method since correction terms were 
not used, and only a relatively few points were used to determine 
these slopes. Since the main purpose in obtaining these values 

was to check the approximate magnitude of those values obtained from 
the distribution functions, great significance is not attached to 
these small differences. Even though absolute magnitudes of the 
densities computed from the distribution functions are not certain, 
the magnitudes relative to each other are thought to be quite good 
and this fact is used later to correlate density with the predicted 


Langmuir frequency and to make a statistical estimate of the 


“The trapezoidal rule gives an approximate value to the area under 
an oa oat ise. ,whwea = 4(y,/2 a ey ee y,,/2) - 
In our integration we used a A of two electron volts giving a fair 
approximation to the integral. Using smaller increments gives an 
increasingly better result. Simpson's rule is more accurate. We used 
the trapezoidal rule because it is suitable for making a running tabu- 
lar integral. Simpson's rule is not conduCive to this treatment. 
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standard deviation for the values of density obtained. 

(e) A plot of relative electron density vs X is shown in 
Fig. 31. These data were obtained from ln 1 Vs . curves by 
using the current at space potential and the electron temperature. 

(f) A plot of relative electron density of high and inter- 
mediate energy beam electrons is shown in Fig. 32. These data 
were obtained by Waymouth's technique of analysis of _ vs Vp 
probe characteristic curves to obtain the energy of the groups. 
Medicus' technique, based upon a diffusion law which is dependent 
upon pressure, was used for relative magnitudes. Our pressure 
range was significantly different (lower) from that of Medicus’ 
development; hence, the results are good only to an order of 
magnitude. 

(g) A plot of eV (average electron energy) vs X is shown in 
Fig. 33. These data were obtained from the average energy of the 


two electron groups. The average energy of each group being 


obtained from the relation 





eV5 
ye’? a°t 
Vv = d(eV 
‘ eV) 2 —) 
dV, 
eV 2. 
V2 qd“¥ 
FP dev) 
eV. a wa 


t » 


These values are shown on the distribution functions in Figs. 19 


through 28. 
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Oy splot of T, (electron temperature) vs X is shown 
in Fig. 35. These data were obtained from the slope of the 
linear portions of the ln iy vs Vp curves after the current 
attributable to the high energy electrons had been subtracted out. 

(i) A plot of T, vs X for both high and low energy electron 
groups is shown in Fig. 36. These data were obtained from the 
average energy of the two electron groups by using the relation 
eV = 3/2 kT, tages is strictly true omly if the disterburciom is 
Maxwellian. The distributions obtained were close to this for the 
low energy groups and the relation was used for the high energy 
groups as well for want of a better relation to use. This plot is 
essentially the same as Fig. 33, since both were obtained from 
the same data. The ordinates differ only by a constant. 

Gr 2 plot of Vp vs X is shown in Fig. 37. These data 
were obtained by taking the maximum of the first derivative for 
positive values of Vo as being space potential. This corresponds 


to the inflection point of the 1 vs Vy curve, or the first devia- 


tion from linearity on the ln i, vs Vp jodKe nus 


p 
(k) A plot of Vp vs X is shown in Fig. 38. These data were 
obtained by taking the first deviation from linearity on the Jn Ip 


vs Vp plot as being space potential. Agreement between this figure 
and Fig. 3/7 is quite good. More faith is placed in the values ob- 
tained from the first derivative because these values are easier to 
read precisely than the deviation from linearity on the logarithmic 
pot . 

(1) A plot of E (electric field strength) vs X is shown in 


Fig. 39. These data were obtained by taking the changes in space 


potential as a function of distance (X) from Fig. 37. 
6 





Plasma oscillations in two frequency ranges were found. 
Kilomegacycle signals detected are believed to be electron 
oscillations and kilocycle signals are believed to be ion 
oscillations. 

High frequency oscillations of 472 Mes and 520 Mcs were 
found at positions 20 mm and 24 mm, respectively, from the cathode 
for the 6.9 micron, 50 ma discharge. Computing the density by use 
of Langmuir's theory relating the density to the slope of i vs 
Vo curves, we obtain the plot of total density shown in Fig. 304, 
below. The densities corresponding to the observed frequency, 
Calculated by using Langmuir's relation f = 8980 nee are shown 


as dashed lines on this plot. 


A es 


520 Mcs 


Electron Density x19? (em73) 





X (mm) 


Pig. 30a. Electron density as a function of cathode - probe 
distance. Data obtained from slopes of i vs V, plot 


p 
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Neither of these densities, corresponding to the observed 
frequencies, correspond with measured points on the total density 
plot, but instead, represent plasma oscillation frequencies lower 
than would be obtained from total density values. 

However, if we make use of the density plot obtained from the 
energy distribution functions and also use Waymouth's technique to 
obtain the distribution of high energy beam electrons, we obtain 
the plot of densities shown in Fig. 29a, below. By using these 
techniques we are able to break our electron concentrations down 
into groups which are: high energy beam electrons (obtained by 
Waymouth's technique) (see p. 23, Section 2, Previous Work), 
intermediate energy beam electrons (obtained from energy distri- 
bution functions) and random plasma electrons (obtained from 


energy distribution functions). 
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Fig. 29 a. ® Total electron density (from distribution functions) 
O Plasma Density (from distribution functions). 
xX Intermediate energy beam electrons (from 
distribution functions) 
— High energy beam electrons (Waymouth's technique) . 
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If we apply the relationship. cited by Kharchenko, et al. 


(31) (described on page 26 , Section 2, Previous Work) . 


(1) f= ¢£ 


O p 
Des 
— 
2 
V 
O 


where S is the thermal plasma velocity and Ve is the beam 


electron velocity and put this in terms of energy, we have 





(2) ae = fy 
L 
~_ 2 
l - eV, 
a 


where eV, is the average thermal plasma energy and eVj, is the 
average energy of the beam electrons. 

We know f,, av and eV, experimentally. Let us now calculate 
the equivalent density corresponding to fy and see where this fits 


em tne density plot, Figs 29a. 


Consider first i = ocoeMes- 


) eV, 


Senye2 =8025 


©Vbintermediate 





ie 
(4) and, if (1 - 0.25)° 520 Mes 


PH 
i 


folehes, 


% 
Realizing this was developed for ie S. 
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The equivalent density is 
oualps 9 , 3 
(S) rn, ={ 451 x 10 Seer It particles/em « 
va “Sau 


Gagilanly;, for f, = 4/2 Mes, 





(6) eV, 
= oc 
eV 
Dhigh 
aM 
Oem i = A- 0.1)° 472 Mes 
i = 447 Mes. 
p 


The equivalent density is 


= 10 particles 

Reterring to Fig. 29a, densities nj) yand m5) (and their 
corresponding frequencies) are indicated by dashed lines. It is 
apparent that these densities correspond with the curve repre- 
senting the plasma density, not total density. 

In obtaining the data for Fig. 29a,we were beset by a 
dilemma. We have previously mentioned that the densities ob- 
tained by the second derivative method were in error by some 
constant factor caused by improper machine calibration. However, 
we are confident that the relative values of these densities are 
rather good with respect to each other. Now the absolute values 
of densities obtained by the Langmuir method are probably accurate 
to about 10%. In order to use our density data obtained by the 
second derivative method to best advantage we have determined a 
scale factor to be applied to these densities (obtained by the 


second derivative method) by taking ratios of each corresponding 
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ordinate of Pig. 29 (second derivative data) and Pedi ree) 
(Langmuir data). We have used the average of these ratios 


to correct both curves to the same absolute magnitude. 


Then, by statistically analyzing the ten density ratios so 
obtained, we are able to compute a number which gives a certain 
measure of confidence in the values of density plotted in Figs. 
29a and 30a. We call this number s to denote a standard devia- 


Pwratio, and plot is as “9° = + s°n_ particles/ em? OT sialas:: 


F 
29a and 30a, for the total density and plasma density curves. 
eewds Calentated to bew+ 0.17. N_ p is used to denote any ordinate 
to the density plots and o the standard deviation in particles/ 
em? . 
Now it is quite interesting to note that neither density 

corresponds to the Langmuir frequency within the confidence 
limits when we attempt to infer it from the total density plot 
However, when we explain our observed frequencies by correlating 
them with the plasma density curve, we note that these densities, 
ny and ng , fall within one standard deviation of this curve, if 
we consider that part of the plasma density Curve at the point of 
disintegration of the beam density Curve. 

The conclusions to be drawn from this analysis are brought 
out in the next section of this report. 

It is worthwhile to note that only the 472 Mcs oscillation had 
a large amplitude (2 ma of detector current). The 520 Mcs oscilla- 


tion was much smaller in amplitude (0.01 to 0.02 milli-micro amp 


of detector current). We associate the strong 4/72 Mcs oscillation 
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with the intermediate energy beam electrons (eV, = 22 electron volts) 
and the’ 520 Mes oscillation with the high energy beam electrons 

(eVp = 50 electron volts). Note also that the density of this 
group of high energy beam electrons is about one order of magnitude 
smaller than either the plasma electrons or the intermediate energy 
beam electrons. 

Other high intensity electron oscillations were detected at 
fener @urrents much more readily than at 50 ma “discharge current. 
Unfortunately, this was discovered so late as to preclude taking 
correlating probe data. The results of these measurements, for 


one probe position, are shown in tabular form below. 


Discharge Current Observed Approximate™ Intensity of 
(ma) Frequency Computed Signal 
(Mcs) Frequency (Detector 
(hie s) Current - am 
100 703 710 3 x 19714 
105 704.5 720 1-2 x 1077 
110 708.5 730 ux 107/ 
115 713 737 ux 10” 
ILO 713.3 746 eee na 
125 Ties 752 cpvouil 


For this probe position no oscillation signals were detected 
below 100 ma discharge current. This suggests that a threshold 
effect for interaction may be present as a function of current 
density. 


* 
By interpolating in Fig. 14 using a pressure correction factor. 
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Fig. 40 is a plot of current density (Jg) vs observed fre- 
quency squared (£5) Showing the linear relation of these values. 

The other type of oscillation found was a low frequency sig- 
nal of about 100 Kes. The ion frequency predicted by calculation 
from the density of the intermediate energy beam electrons is 
135 Kes. Hence, these oscillations are believed to be ionic 
in nature. They existed throughout the length of the tube (for 
the 50 ma discharge at 6.9 microns pressure), maintaining 
essentially the same frequency; however, the intensity of the 
Signal did vary with position in the tube. The results of these 
intensity variations are shown in Fig. 34, where I (intensity in 
volts) is plotted vs X. It is also significant to note that the 
oscilloscopic traces of these oscillations show a more fixed 
adherence to a steady fundamental frequency and waveform on the 
cathode side of the discharge where the high and intermediate 
energy electron beams, coming from the cathode, were in evidence. 
On the anode side of the discharge (in the region after which 
scattering has occurred) these oscillations show a wider spectrum 
of frequencies (varying between 90 and 110 Kes) and the intensities 
of signals and waveshapes are quiie variable when sampled in time. 
Fig. 41 is an oscillographic photograph of an oscillation pattern 
on the cathode side of the meniscus region. Fig. 42 is an oscillo- 
graphic photograph of an oscillation pattern obtained on the 


anode side of the «wniscus region. 
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Visual and photographic observation of the discharge was 
made in effort to try to discover a meniscus, a scattering region 
and/or a diverging beam of electrons. No meniscus was observed 
at low currents; but, visual observation seemed to indicate a 
region of maximum disturbance in the beam, and at about the same 
location, a tendency for the beam to diverge from the shape dis- 
played in emerging from the cathode. The region in which this 
phenomenon appeared is indicated on Figs. 29 through 39 as the 
eross-hatched region. Figs. 43 through 45 are photographs and 
sketches of the discharge which show the beam, its tendency to 
diverge, and apparent meniscus formation. Note that meniscus 
Formation appears to grow with increasing current. 

In attempting to correlate frequency measurements made at 
higher currents during the second data run (50 ma discharge at 
6.9 microns of pressure) with the current-density data of the 
first run (five different currents at 5.9 microns of pressure) , 
an interesting fact (overlooked as being trivial when it occurred) 
was noted. Referring to Fig. 14 for a discharge current of 125 ma, 
at a position ten mm from the cathode, a peak density occurs. The 
minimum of this gradient occurs at 12 mm from the cathode. Referring 
to the original probe characteristic curve corresponding to this 
point (Fig. 7, 125 ma discharge current), it was noted that consid- 
erable irregularity of the ion portion of this curve was observed 
in a probe voltage range of -40 to -100 volts with respect to the 
anode. The phenomenon manifested itself in a step increase in 


probe current of about 0.004 ma which corresponded to a step 
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increase of about 0.4 volts. This distortion in probe charucter- 
istic lasted out to about - 100 volts at which point it tended to 
return to its normally expected shape, but with some noticeable 
change in slope. At the time of occurrence, the peculiarity was 
noticed and the curve retraced several times, each time giving a 
similar result, although the exact voltage at which the step in- 
crease in current occurred did not coincide exactly, it held to 
within a 20 volt region. Sweep speed of probe voltage was in- 
creased with no change in effect -- in fact the speed of sweep 
was increased to about 50 cps (vs the normal 0.05 cps) and the 
trace photographed on the oscilloscope. The same sort of 
phenomenon showed thereon except that this trace seemed to show 
oscillating signals riding on the probe characteristic from about 
-20 volts out to the maximum excursion of -130 volts. Further 
examination of the electron portion of the probe characteristic 
curve disclosed that this was the first probe characteristic curve 
(for this current), in moving the probe from cathode to anode, in 
which a section having a constant slope in the positive current 
region (denoting a region separating two energy groups, ref. 
discussion of Waymouth's work, p&. 23 , Section 2, Previous 
Work) did not eae 

A similar occurrence took place for the same discharge current 
at 10 mm from the cathode. This disturbance to the probe was less 
pronounced, amounting to a change in slope on the linear 1p, vs V 


p 


Curve and a distinct deviation from linearity on the Cae Vs Vp curve. 


“See Fic. 46 whiech® 1@lustrates this. ci Fig. 85 in which only 
one electron group is present. 
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Still considering this 125 ma discharge, it was noted at 
the time of taking the data that the position 10 mm from thi 
cathode was just barely on the cathode side of the diverging 
beam region (observed visually) of the discharge. The electron 
portion of this probe characteristic curve showed a constant . 
slope region denoting two distinct groups of electrons still 
present at this position in the tube. The position 12 mm from 
the cathode was inside, i.e., on the anode side, of the diverg- 
ing beam. The meniscus is believed to have been present about 
10 mm from the cathode. 

Similar occurrences were evident on the probe character- 
istic curves at 14 mm from the cathode for the 100 ma discharge and 
at 22 min (see Fig.8) fror tte cathode for the 75 ma and 125 ma, dis- 
charges at probe voltages of about -100 volts with respect to 
the anode. These points corresponded with a peak density 
gradient in each case (see Fig. 14). 

In light of what was discovered on the second data run 
(50 ma discharge at 6.9 microns pressure) it now appears clear 
what happened in these cases during the first data run. These 
disturbances to the probe characteristic are associated with the 
same region in which (a) electron oscillations were detected, 

(see Fig. 29). (b) ion oscillation intensity reached a maximum 
and ion oscillation waveshape became most distorted, (see Tig. 34). 
(c) beam electron concentration largely disappeared and merged into 
one random distribution (see Figs. 23 and 24), and (d) either a 


maximum or a minimum density gradient was present {see Fig. 14) 
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The changes in slope as far negative as -100 volts ace 
believed to be caused by electrons acquiring enough energy in 
the interacting region in the vicinity of the meniscus (where 
oscillations are believed to originate) so that they can overcome 
this amount of retarding field and become collected by the probe 
There are, it appears, distributions of electrons in these inter- 
acting regions having maximum energies as high as 100+ electron 
volts. The step increase in current shown on Tig. 7 is believed 
to be a monoenergetic group of electrons having an energy of 
about 50 electron volts with a variability of about + 10 electron 
volts when sampled in time. 

A distinction can be made with respect to the effect on the 
probe characteristic as it is affected by the density gradient. 
A large disturbance, i.e., a step increase in probe current, and 
a small disturbance, i.e., a change in slopeywas noted where the 
position was associated with a minimum of the density gradient 
(corresponding with Emeleus' findings that oscillations reach a 
maximum intensity at the bottom of a density gradient}. Only 4 
small disturbance, i.e., a change in slope of the probe character- 
istic curve, was noted where the position corresponded to a muxi- 
mum of a density gradient (in agreement with Emeleus' findings 
that oscillations are initiated at or near peaks in the density 


gradient). 
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The effect of probe contamination was studied to determine 
whether cathode sputtering was causing the probe to coflegt im- 
pure material and thus change its collecting properties. A 
curve was traced with a freshly cleaned probe (by ohmic heating 
to a cherry-red condition) then allowed to sit immediately in 
front of the cathode, perhaps 2 mm away, (the worse possible 
position), for a period of 15 minutes to one half hour and then 
the same characteristic curve was retraced. A siight change in 
the characteristic was noticeable; however, when the interval of 
time was reduced to the same time as was used for making a set of 
data curves, the effect of probe contamination on the character- 
istic curves was entirely negligible. This is probably attribut- 
able to the very low evaporation rate of material from the cathode, 
given as less than 0.8 micrograms of barium per em* per hour by 
the manufacturer. 

Miscellaneous tests were conducted to determine the effects 
of electronic errors associated with the measurement of probe 
characteristic curves. Probably the best estimate of the reiia 
bility of the probe characteristic curves was obtained by aillou- 
ing the curve to be retraced 15 to 20 times to note the driit 
tendency due to the amplifiers and associated electronic circuitry 
Results of these measurements showed that the maximum deviation ta 
be expected was of the order of + 0.002 ma of current for electron 
currents of the order of g.l ma. Curves were also re-run after an 
elapsed time of days and no significant change in characteristic 
was noticed except when the gas was obviously dirty because of 


long periods of operation (up to several weeks). 
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6. Conclusions. 

The results of our second data run (50 ma discharge at 6.9 
microns pressure in argon) as shown in Fig. 29 tend to substanti- 
ate that plasma electron oscillations originate at the peak of a 
density gradient and grow in amplitude down a density gradient as 
originally proposed by Allis (30) and observed experimentally by 
Emeleus and co-workers (33). 

Electron oscillations were detected only at positions in the 
tube immediately after divergence of the electron beam had oceurred. 
These oscillations could be associated with a peak density gradient, 
and the observed frequency of these electron oscillations bore a re- 
lationship to the predicted Langmuir frequency when correlated with 
the plasma density (vs' the usual total density). The reason for 
this will become clear later in discussing a possible mechanism 
mor oscillation excitation. 

Ion oscillations were detected throughout the length of the 
tube with a fundamental frequency of about 100 Kes. Fig. 34 shows 
that the intensity of these oscillations fell to a minimum just 
prior to the point where the total density gradient reached a 
maximum (see Figs. 29 and 34) and then the oscillation intensity 
rose rapidly to a maximum where the total density gradient was a 
maximum. This region is presumed to be that where particle inter- 
action is greatest and where energy exchange is taking place and 
thus exciting electron oscillations. This region is outlined by 
cross-hatching on Figs. 29 through 39. The fact that the photo- 


graphic traces of the ion oscillations across this region, in 
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going from cathode to anode, show considerable regulurity on the 
cathode side (see Fig. 41) and more intensity and more irregularity 
in wave shape and frequency on the anode side (see Fig. 42) 
indicates that in this region, strong microfields are present 

which are disturbing the ion oscillations. The change in intensity 
of the ion oscillations indicates that ions are taking part in the 
energy exchange mechanism by losing energy as they move out of the 
turbulent region toward the cathode. 

Disturbances to the probe characteristic curves noted on the 
first data run being associated with maximum or minimum density 
gradients points strongly to the conclusion that strong fields 
and oscillations were present, and considerable energy exchange 
was taking place, at these positions in the tube. Disturbances 
out to -100 volts suggests that some electrons were acquiring 
better than 100 electron volts of energy These fast electrons 
must be generated in the interacting region since the maximum 
energy that could be gained from the discharge tube would be 
of the order of the tube voltage drop, about 20 electron volts. 

Maximum disturbances were noted at positions corresponding 
to a minimum of the density gradient in agreement with Allis‘ 
contention that oscillations grow in intensity down a decressing 
density gradient. Smaller disturbances to the probe character- 
istic curve were noted at positions corresponding to a maximum 
of the density gradient in accord with Allis’ idea that plasma 
oscillations originate at, or near, a peak density gradient. At 
higher discharge currents, multiple density gradients, (i.t&., two 


peaks) were present (see Fig. 14). Probe disturbances were 
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associated with each of these. This fact suggests that a more 
complex interaction is taking place with perhips the formition 
of two such interacting regions, one following the other 

Examination of the distribution functions shown in Figs. 19 
through 28 show that two groups of electrons are present on the 
cathode side of the interacting (meniscus) region. On the anode 
side of the interacting region, only one group of electrons is 
present. Further, by employing Waymouth's analysis of the probe 
characteristic curve, a third higher energy group of electrons is 
also found to be present, but these electrons are either much less 
numerous or of too high an energy to show on the energy distribu- 
tion function. We have earlier denoted these higher energy groups 
as an intermediate energy group of beam electrons and a high energy 
group of beam electrons, respectively. 

These beam electrons emerge from the region in front of the 
cathode and interact with the plasma electrons at a region in the 
tube where visual divergence of the beam and the aforementioned 
disturbances in the probe characteristic curves occur. 

The evidence of er data runs has suggested to us that the 
mechanism exciting the oscillations is associated with these beam 
electrons, because all of these associated phenomena are indicutive 
of a region of instability in which local fluctuating microfields 
of high strength would be likely to exist. 

In the light of this evidence, the authors believe that the 
beam electrons, upon reaching the plasma, are strongly interacting 
with the plasma (when beam current densities are adequate), pro- 


ducing local instabilities which excite plasma osciltations. 
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Ion oscillations may well play a role in the excitation process, 
or may only be a by-product of the electron beam excitation 
mechanism. 

The theory of such a mechanism has been developed by Akhiezer 
and Fainberg and is discussed by Kharchenko, Fainberg, Nikolayev, 
et al., (31). The important points in their discussion are sum- 
marized in Section 2, Previous Work, page 26. 

As previously mentioned, there are two possible types of 
instabilities, convective and alsolute. In our discharge, con- 
vective instability is believed to have occurred since oscillations 
were detected only in a very small region associated with the 
initiation of the instability, and none were detected on the anode 
side of the point of initial detection where one would expect the 
oscillation intensity to increase as the oscillation moved down a 
density gradient. Possibly the spread in velocities as the beam 
electrons became randomized in their initial interaction with the 
plasma tends to remove the instability as suggested by Kharchenko, 
et al., (31). This may be the reason that oscillations are not 
detected other than in the region of initial scattering. 

Perhaps as a bit of speculation, we might mention that the 
possibility of modified "runaway" electrons (31, 1, 36) being 
generated and initiating the instability should not be entirely 
ignored. Although “runaway™ electrons have not been observed 
(to our knowledge) in discharge tubes such as ours, conditions 
suitable for their production appear to be met in our tube and 
some of the phenomena observed might well be attributed to their 


presence. For example, the sudden termination of the instability 
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and damping out of oscillations are what one would expect if the 
excitation was produced as a result of "runaway" electrons. 

Or possibly, electrons traversing the internucting region may 
be given a pulse of energy by the oscillations and emerge as 
"runaway" electrons. These “runaways" would not interact with 
other particles because their collision prolbwéhity is mech 
less at their greater energy, thus, they would completely emerge 
from the plasma and finally give up some of their excess energy 
by collision with walls of the tube. 

One other point is worthy of discussion at thi» tim:. On the 
second data run, oscillations were detected at higher dischurge 
currents having signals of considerably greatéry amplitude than 
those observed at 50 ma discharge current. All efforts at detecting 
oscillations at discharge currents lower than 50 ma were unsuccees-. 
ful. Examination of the density plot from the first data run 
(Fig. 14) shows that no sharp gradients exist for the 25 and 50 
ma discharges. Sharp gradients do exist for 75 ma and above. In 
the light of this observation and the experience obtained on the 
second data run, it is believed that current density (in order to 
produce a sufficiently strong interacting beam), and pressure are 


important in establishing the instability and hence osciliutions. 


i 
PI 


(3) 


In particular, had we operated at slightly higher pressures 
and/or slightly higher currents, it is believed that oscillations 
would have been found more readily and that a closer approach to 

an ideal plasma (i.e., having a Maxwellian distribution of energies 
would have been obtained since a smaller mean free path would have 


produced more randomization of electron energies. 
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The possibility also exists that we were operiting in a 
current density - pressure region which was jueét adequate to allew 
a beam to form and move some distance before interacting with the 
plasma. Operation at higher current densities has shown us that 
the meniscus region moves toward the cathode as «umrrent if increased 
Further experimentation should reveal the optimum currert censity 
pressure relationship to use in attempting to isolate the excitation 
mechanism. 

The results of this investigation are based upon only Limited 
data, some of which is questionable in magnitude (perhaps: 10 to 15 
per cent in error in ordinates of the energy distribution functions} 
because of the equipment troubles at the end of our work. Howe vex, 
all data were qualitatively good, at least. Much more dit.u need: 
to be collected before definite, positive conciusions may be drawn. 
Nevertheless, the correlation between all of the evidence pre sented 
is quite convincing and seems to substantiate the tentative con- 


clusions which have been drawn. 


cs 2 
The area of the cathode used in this work is WT cm*; hence, 
dividing any of our discharge currents by fT will give the current 

density. 
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7. Recommendations. 


In keeping with our findings and observations, the foliowing 


recommendations for further work are made: 


(a) 


(b) 


(c) 


(d) 


Continue probe studies and oscillation detection ar 
higher pressures and/or higher currents, perhaps 
varying pressure aS a parameter to determine chang#s 
in the form of distribution functions. 

Establish the pressure~current density range in which 
oscillations are detected most readily, then attempt 
to find a simple oscillation-density gradient rela- 
tionship which will present the least complicated 
Strueture initialing 

Once suitable oscillations are found, make probe 
studies and frequency measurements of both ion and 
electron oscillations to determine where in the tube 
they arise and attempt to correlate these measure- 
ments with visual observation, energy distribution 
functions, and perhaps also some spectroscopic 
observations. 

Study the region is which osciilations appear to 
arise in great detail; i.e., make probe and oscilla- 
tion measurements around this region utilizing very 
small increments of movement of the probe, say 0.5 
mm or less, in an effort to isolate the mechanism 


of excitation. 
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Figure 3 


Power Supplies and Probe (Cj 
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Figure 29. Electron density as a function of Cathode - Probe rosition. 
Obtained from energy distribution functions. 

® -- Total density 

O -- Density of plasma (low energy group) 

MX-- Density of beam electrons (intermediate energy group) 
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Figure 30. Total electron density as a,function of Cathode - Proke 
position, Obtained from the slope of 1a vs. V_ curves. 
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Figure 31. Relative total electron density as a function of Cathode - 
Probe position. Obtained from In - VS. "5 CULVES neat Space potertial. 


-“-~ 
a 

F Argo 

~ eL/ 2 6.9 microns of mercury 
a 

LW 

a o 

Vv) 

2 2 

a) 

a 

& 

c 

O & 0 

MH 

aed 

5 | 

wo | 

rd 

ra) 

> 

ord 

LW 

o ; ) 

v 

; AN 

0 8 16 L, oS 40 
C A 


X (mm) 


Figure 32. Relative electron density of beam electrons 4s a function 
of Cathode - Probe distance. 

O-- Intermediate energy group 

X-- High energy group 
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Figure 33. Average energy of electrons in electron volts as a functior 

of Cathode - Probe distance. Obtained from energy distribution functions. 
-- Intermediate energy group . 
-- Low energy group 
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Figure 34. Ion oscillation intensity as a function of Cathode = Probe 
_ distance. 
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Figure 35. Electron temperature in degrees Kelvin as a function of 
Cathode - Probe distance. Obtained from ln 7 VS. : cutves. 


A -- High energy group In i vs. V_curves corrected for ion current 
@ -- Low energy group and eee Rnto a high and low energy group. 


O-- Low energy group ) Uncorrected linear portions of i vs. V_ curves 
used for these data. P P 
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Figure 36. Average electron temperature in degrees Kelvin as a function 
of Cathode - Probe distance. Obtained from energy distribution functions. 
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Figure 37. Space potential in volts as a function of Cathode - Probe 
distance. Obtained from maximum of first derivative of Lo V8.0 CUIVES. 
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Figure 38. Space potential in volts as a function‘of Cathode - Frote 
distance. Obtained from the first deviation from linearity of Im i 
vs. V5 curves. P 
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Figure 39. Electric field strength as a function of Cathode = Probe 
distance. Obtained from data of Figure 37. 
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Figure 40. Current density as a function of the square of the observed 
requency of electron oscillation for the probe located immediately 


after the region of suspected maximum interaction. 
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Fig 4la. Several oscilloscope traces 
of ion oscillations. 








Fig 41b. Oscilloscope trace of ion 
oscillation when sampled in time. 





Fig 4lc. Oscilloscope trace of ion 
oscillation when sampled in time. 


Fig 41. Oscilloscope traces of ion oscillations on the 
cathode side of the meniscus region. Discharge current 

50 ma. Pressure 6.9 microns of mercury. Cathode-anode 
sracing 40 mm. Probe position 12 mm from cathode - 
oscilloscope sweep speed Sygsec/cm. Amplitude 1 volt/cm. 
Wote smaller amplitude, more fixed waveform and more steady 
frequency compared with Fig 42. 
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Fig 42a. Several oscilloscope traces 
of ion oscillations. 
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Fig 42b. Oscilloscope trace when 
sampled in time. 





Fig 42c. Oscilloscope trace when 
sampled in time. 


Fig 42. Oscilloscope traces of ion oscillations 

on the anode side of meniscus region. Cathode-anode 
spacing 40 mm. Probe position 20 mm from cathode. 
Discharge current 50 ma. Pressure 6.9 microns of 
mercury. Oscilloscope sweep speed 5 ysec/cm. 
Amplitude 1 volt/cm. Note the greater amplitude, 
distortion of wave form and spreading of fundamental 
frequency compared with Fig 41. 


99 


aa) OO a be eee) aE gre 
bay’ 70s ey 


Warn ee quent eee wes 2 )¢ 
ae A Mm OS” ieee 








Figure 43 a. Photograph of 50 ma discharge 
in argon at 6.9 microns of mercury. 
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Figure 43 b. Sketch of 50 ma discharge in 
argon at 6.9 microns of mercury. 
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Figure 44 a. Photograph of 120 ma discharge 
in argon at 6.9 microns of mercury 
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Figure 44 b. Sketch of 120 ma discharge in 
argon at 6.9 microns of mercury. 
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Figure 45 a. Photograph of greater than 150 ma 
discharge in argon at 6.9 microns of mercury. 
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Figure 45 b. Sketch of greater than 150 ma 
discharge in argon at 6.9 microns of mercury. 
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Fig 47a. Energy Distribution Function. Argon 
. Current = 50 ma 
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N, = f NG(eV) d(ev) => | NP(ev), (AeV),, « If Lev, = eV, Lev 
and we use the trapezoidal approximation to the area, then 
oo =16 d= 28 
Ny 2 Aov| ang (ev) + > NB (ev), N, 4 Aev| > | NB (eV), + SNZ (eV), 
oC = X<216 
_ MEho a : 
N. £ 1.722 x 1023 (Aev)| ava" ve it, | , where i = d*i 
4 eVxa2 Pp 
"as dv 
Since we are using 1/3 i", we multiply by 3 and factor p 
out 107° from V“i" terms and our compuccetgoe: formula becomes; 
7 K=16 6 
Ny = 1.722 (3)(2 eV) x 10°] yo/2 + > y,|, where y = ve! x 10 
kan ~ 


8 &=16 
Ny = 1.033 x 10 | y2/2 + 7h ve 
x=2 


It is convenient to express this computation in tabular form 
in which a running integral may be set up and the integral 
eV 
2y N@(eV) d(eV) computed at the same time for the purpose 
eV 
of fevaluating the average energy, eV, of a group of electrons, 


Figure 47. Example of electron density computation and average 
energy computation from energy distribution function. 
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16 6 
z fv N@(eV) d(eV) )* i’ d(ev) 
1.@ey ev, fea 0 = K - 
; . Pa 
fxpvev) d(ev) a i" d(eVv) 
By the same process, No and eV. may be computed. 


30 
Then, Ni, = Ny + Ny = J MW) d(ev) 


These computations are shown in a tabulation on the 
next page in Figure 47b. 


Figure 47. Example of electron density computation and average 
energy computation from energy distribution function. 
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